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Abstract 



We present a new simulation engine for fast generation of ultra-high energy cosmic ray data based on parameterizations 
of common assumptions of UHECR origin and propagation. Implemented are deflections in unstructured turbulent 
extragalactic fields, energy losses for protons due to photo-pion production and electron-pair production, as well as 
effects from the expansion of the universe. Additionally, a simple model to estimate propagation effects from iron 



nuclei is included. Deflections in galactic magnetic fields are included using a matrix approach with precalculated lenses 
generated from backtracked cosmic rays. The PARSEC program is based on object oriented programming paradigms 
enabling users to extend the implemented models and is steerable with a graphical user interface. 
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1. Introduction 

Recent results of the Pierre Auger Observatory imply 
that the ultra-high energy cosmic rays (UHECR) are ac- 
celerated at extragalactic sources with a distribution fol- 
lowing the large scale structure [HE]. Nevertheless, the 
exact origin of these cosmic rays remains so far unknown. 
UHECR are likely charged particles [2j [3] and thus de- 
flected in extragalactic and galactic magnetic fields. They 
can therefore be considered messengers of cosmic acceler- 
ators as well as fields they are subjected to while propa- 
gating to the Earth. During their propagation from the 
sources to the Earth UHECR lose energy due to inter- 
actions with different photon backgrounds. In the case of 
protons the dominant processes are photo-pion and electron- 
pair production, while for nuclei also photo-disintegration 
has to be accounted for. The composition of UHECR is 
still under debate [4, 5 . 

Progress in the understanding of origin and propaga- 
tion of UHECR is enabled by comparison of measured 
data with simulated UHECR data. A second field of ap- 
plication for simulated UHECR data is development and 
benchmarking of observables suited for disentangling spe- 
cific aspects from the observed data. In this publication we 
present a Monte Carlo generator for UHECR propagation 
from extragalactic sources to the observer. 

As the charges of the UHECR, their origin, and the 
magnitude of their deflection in cosmic magnetic fields are 
linked, interpretation of cosmic ray data relies on the com- 
bination of four different types of models: 
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1. A model describing position, spectrum, and relative 
luminosity of the sources, 

2. A model for the deflection in extragalactic fields, 

3. A model describing the composition and thus the 
corresponding energy loss of the particles, and 

4. A model describing the deflection in the galactic 
magnetic field. 

To test these models a fast generation of sufficiently 
large data samples at every point of the parameter-space 
is required. 

The most obvious approach for a UHECR Monte-Carlo 
generator allowing detailed simulations of individual phys- 
ics processes is to follow the trajectories of individual par- 
ticles from the source to the observer and account for en- 
ergy losses and deflections during the propagation. This 
forward-propagation is for example implemented in the 
CRPropa program [6 . In this approach the, compared to 
intergalactic distances, small size of the observer makes the 
generation of large Monte-Carlo data sets, needed for ex- 
tensive parameter scans, highly challenging. As the chal- 
lenge becomes more and more complex for contributions 
from distant sources a first approach combining a full sim- 
ulation with the PARSEC program has been presented 
in 0. 

A different approach aiming at UHECR mass produc- 
tion is to backtrack particles starting at the observer and 
associate them to the objects in the source model. This 
approach is well understood and documented in the lit- 
erature (e.g. [8j|9]). However, it requires huge trajectory 
databases (e.g. [10]) and the trajectories are associated to 
sources only at the end of the simulation. 

To enable UHECR simulations for arbitrary source mod- 
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els with sufficient speed for extensive parameter scans we 
developed PARSEC as a parametrized simulation engine. 
PARSEC calculates the density function of the probabil- 
ity distribution to observe a cosmic ray in a certain energy 
range coming from a discrete direction (pixel) based on 
a distribution of individual point sources. For the calcu- 
lation it uses parametrizations for the energy losses and 
the effects from deflections in extragalactic fields. Pre- 
calculated models for the galactic magnetic field can then 
be applied in a separate simulation step using matrix tech- 
niques. This allows for independent combination of models 
for extragalactic and galactic propagation without redun- 
dant calculations. From the resulting probability densities, 
simulated data sets can be generated quickly. 

This publication is structured as follows: The imple- 
mented models and details of the implementation are pre- 
sented in section [2j The resulting energy spectra, parti- 
cle horizons, mean deflections, and exemplary probabil- 
ity density maps are presented together with performance 
benchmarks in section 03 The conclusions in section [4] 
are followed by an appendix explaining uncertainties in 
the generation of the galactic lenses and optional improve- 
ments concerning the computing time. 

2. Methods 

We interpret the probability to observe a particle from 
a discrete direction ((/>, 6) j on a pixellated sky as an ele- 
ment j of a vector p ohs resembling the probability distribu- 
tion. For every energy range E\ < E^ < E\ a separate vec- 
tor is calculated, leaving the total probability distribution 
V from which individual cosmic rays are later generated 
as a set of vectors V Q bs = {Po&si- 

For every energy range i the calculation of the ob- 
served probability vector p % ohs is performed in two steps: 
First, probability vectors p l for the expected distribu- 
tion including all effects except the galactic magnetic field 
are calculated. In a second step this probability vector is 
transformed by a matrix derived for the galactic magnetic 
field model, yielding p z obs . 

2.1. Extragalactic Propagation 

The discrete probability distribution p l is calculated 
by the sum of the contributions of every individual source 
Sk of a source model {S\ • • • Sn} to every pixel Pj. We 
separate the contribution of Sk to p* into three factors: 

1. A factor fs including the source distance and lumi- 
nosity, 

2. A factor describing the source energy spectra and 
energy loss effects, and 

3. A factor fs distributing the flux from one source 
on multiple pixels to account for the deflection in 
magnetic fields. 



With these ingredients the probability to observe a particle 
with energy Ei in pixel j can be written as 

P]=^i-J2fs(S k )f B fE (1) 

k 

where Ti denotes a normalization factor ensuring J2i j Pj = 
1. 

2.1.1. Energy Losses 

A particle emitted by source Sk with the energy E in j 
propagated a distance r within the extragalactic magnetic 
field model, and is then observed with energy Ei. The 
probability to observe the particle therefore depends on 
the source spectra, the energy loss of the particle and the 
propagation distance, which is summarized in For 
source spectra following a power law described with spec- 
tral index 7 this corresponds to 

fE = lJ rZg { E ir^j]r - E injj) ( 2 ) 

where the scaling factor (1 + Zg)~ x of the universe at the 
cosmological epoch of particle injection accounts for cos- 
mological time dilation. The range of injection energies 
Einj, 1 /r is calculated with the continuous energy loss ap- 
proximation by numerically integrating 

with the total energy loss length defined as L(E^z)~ 1 = 
Lad(z)- 1 + L vh {E,z)- x . Here L ad (z) = c/H(z) is the 
energy loss length for the adiabatic energy loss by the ex- 
pansion of the universe with Hubble parameter H(z) = 
Ho^y^L m ^o(l + z) 3 + &a,o at the cosmological epoch z. Here 
we use for the current matter density f2 mj o = 0.3 and the 
current dark energy density Qa,o = 0.7 [TT] . 

For the energy loss lengths for interactions in back- 
ground photon fields L p h(E, z = 0) as published in [12J[T3] 
are implemented. As an extension, we also implemented 
a simplistic model for the observed UHECR being iron 
nuclei. For this we also use a continuous energy loss ap- 
proximation with an attenuation length given by the max- 
imum of the nuclei calculated in reference [14]. Here the 
iron nuclei do not disintegrate but keep a chargenumber 
of Z = 26. The propagation of secondary particles is not 
included. As the cosmic rays have maximum range and 
deflections, this model yields a maximum isotropy. From 
the energy loss length of photon interactions at redshift 
z = the energy loss length at z is derived from scaling 
using 

L ph (E, z) = (l + z)~ 3 L ph ((1 + z)E, z = 0) (4) 

to account for the increase of the energy and density of 
the CMB background photons. 
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2.1.2. Scattering around Sources 

The effects of the extragalactic magnetic field are pa- 
rametrized assuming a turbulent field in which the par- 
ticles perform a random walk. In a turbulent magnetic 
field, the flux of a single source Sk is distributed over sev- 
eral pixels pj. If the particles perform a random walk, 
the distribution of the angles aj^ between the direction of 
source Sk and the center of pixel p % - follows a Fisher distri- 
bution [15] . which can be regarded as normal distribution 
on a sphere. The second factor fs of eq. [IJthus reads 



d (k-COS CXj,k) 



(5) 



47r sinh (k) 

with the concentration parameter n. 

The concentration parameter k is related to the root 
mean square (rms) a of the deflection for small angles by 
k = I /a 2 . For small angles the root mean square of the 
deflection of UHECR with charge Z in units of the elec- 
tron charge e in turbulent fields with Dk ^> A can be 
parametrized as 




lEeV J 



(6) 



Here Dk is the current proper distance of Sk to the ob- 
server, B is the strength of the magnetic field, and 



7r / ~~ dk o . 7 , 

A = ^2 / ~ B ( k ) 



(7) 



correlation length of the magnetic field [T6l ? ]. 

This parametrization for the rms of the deflection angle 
does not include energy losses. To derive a first order 
approximation including energy losses, we first derive eq.[6] 
with respect to the source distance x 
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Assuming (dE/dx) • (x/E) being small for the ultra-high 
energies considered here the second term of eq. [8] can be 
neglected. Integrating eq. [8] to the source distance 
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then yields the rms of the deflection for particles from 
source Sk- 

2.1.3. Elongation of Propagation Time 

Due to deflection in magnetic fields the length of the 
trajectory of the particles r = Dk + r from sources in dis- 
tance Dk is elongated by an extra distance r. Neglecting 
energy losses, the additional propagation distance r origi- 
nating from small deflections in the extragalactic magnetic 
field can be written as 
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Table 1: Results of the fit of eq. [13] to UHECR densities in forward 
simulations. 



Parameter 


Value 


Unit 


Pi 


0.1 


nG" 2 EeV 2 Mpc" 3 / 2 


s 


5 




x t 


70 


nG EeV" 1 Mpc 3 / 2 


P2 


9 


nG" 1 EeV Mpc" 1 / 2 


Xi 


140 


nG EeV" 1 Mpc 3 / 2 



with definitions as in eq. [6] [16 . To account for energy 



losses in the extended propagation length we write eq. [10 
for infinitely small dx 



dr oc ( 
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(ii) 



with < x < Dk. Assuming again (dE/dx) • (x/E) to be 
small, the second term of eq. [TT] can be neglected. The 
result can be written as a Riemann sum 



r oc 



N 



B 2 AZ 2 



\( X i 



-i) (12) 



with xn = Dk- Eq. [12] is solved iteratively for every indi- 
vidual source Sk- 

2.1.4- Increase of Particle Density 

The factor fs(Sk) accounts for the relative individual 
luminosity Lk of source Sk and the density of particles 
at the position of the observer in current proper distance 
Dk- If the particles propagate on a straight line, the flux 
from a source is distributed on a sphere with radius of 
the current proper distance Dk of the source resulting in 
fs(S k ) = L k /Dl 

Nevertheless, in the presence of magnetic fields the 
density of UHECR is higher compared with linear prop- 
agation. In some distance to the source all information 
about the origin of the UHECR is lost, and instead of a 
directed random walk the density is described according to 
a Wiener process. To model this transition we simulated 
the trajectories of individual UHECR from one source in 
a turbulent magnetic field with a modified version of the 
CRPropa software [6]. Using g = ^^f~A the density of 
UHECR in distance Dk can be approximatly described by 



fs = 



with 



{l+pi-g 2 ■D k )(l-T)+T-p 2 -g-e 



(13) 
(14) 



and parameters as given in table[l] For 0.1 nGEeV 
< g < 2.5nGEeV _1 Mpc 1/2 the deviation of eq. |13| from 



Mpc 



1/2 



the simulated values is typically below 20%. For a ho- 
mogeneous distribution of sources, the resulting spectrum 
agree also within 20% with the universal spectrum [T7] . 
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2.2. Galactic Fields 

To model particle propagation in the galactic magnetic 
field we neglect energy losses during the relatively short 
galactic propagation. The effects of the galactic magnetic 
field can thus be addressed as magnetic lensing [18J[T9]. As 
there is no random process in this model, a particle with 
energy Ei entering the galaxy at a point C on the surface 
of the galaxy at a discrete direction denoted with index n 
is always observed at a direction denoted with index m. 

The sources considered here are generally in a large 
distance compared to the size of the galaxy, which reduces 
the galaxy to a point in view of the source. The directions 
of entry n for particles with energies Ei can therefore be 
averaged over all points of entry C. Thus a particle with 
energy Ei entering the galaxy from direction n can be de- 
flected into several observed directions m. The probability 
of observing a particle on Earth from direction m which 
entered the galaxy from direction n is / m , n . The / m?n form 
a matrix Ci which represents the galactic lens for energy 
E^ The magnetic lensing can thus be written as a matrix- 
vector equation which transforms the extragalactic prob- 
ability vector pl g as described in the previous subsection 
reading 

£i'P l eg =P l obsr (15) 

The model for the galactic field is hence completely 
described by a set of matrices {C\ • • • jCn} with the energy 
index i = 1 • • • N. 

The individual matrices Ci can be generated by back- 
tracking cosmic rays with isotropic starting directions from 
the earth with the following technique. The starting di- 
rections of backtracked particles are binned in N pixels 
indexed by m. The directions in which the cosmic rays 
leave the galaxy are binned into N pixels indexed by n. 
Counting all trajectories leads to a matrix Ci with ele- 
ments Im.n- We normalize Ci by the maximum of unity 
norms ||£$||i of all lenses reading 



1 



max \ \Cj\ 



(16) 



Each element l mjn of Ci is the probability that a parti- 
cle entering the galaxy in pixel n is observed in direction 
m. This normalization ensures that an isotropic distribu- 
tion outside the galaxy is mapped to an isotropic observed 
distribution as required from Liouville's theorem (cf. [18]). 

In general the galactic field modifies the energy spec- 
trum of cosmic rays depending on the positions of the 
sources Sk as the flux from individual regions in the sky is 
suppressed or enhanced [T9] , While the relative deforma- 
tion of the energy spectrum is accounted for in the normal- 
ization procedure described above, no information about 
a suppression of the total flux by the galactic magnetic 
field is obtained by the generation of the galactic lenses 
from backtracking. Such information could be obtained 
by generation of the lenses from forward tracking or for 
analytically derived lenses. 



Galactic lenses, generated from backtracking Monte- 
Carlo data in the described way, introduce an uncertainty 
in the observed probability distribution. This uncertainty 
is discussed in | Appendix A.l| 

2.3. Technical Realization 

PARSEC is implemented as C++ code with a Python 
interface. It is based on the Physics Extension Library 
(PXL) [20]. PXL is a collection of C++ libraries with a 
Python interface providing classes and templates for ex- 
periment independent high-level physics analysis. The us- 
age of the PXL libraries facilitates modular object-oriented 
programming and allows graphical steering of the simula- 
tion components using the VISPA program [2Tj . 

The individual simulation steps are implemented as 
separate PXL modules which can be individually connected 
and configured to a simulation chain using the graphical 
user interface (GUI) of VISPA. A realization of a UHECR 
scenario is represented by a data container, which is con- 
secutively processed by the following modules. 

Source Model 

Sources of UHECR are represented as individual ob- 
jects. They are added to the realization with user-defined 
coordinates in a Python or C++ module. An exemplary 
module for isotropic source distributions is included in 
PARSEC. Modules generating sources e.g. from astronom- 
ical catalogues can be created by users. 

Extragalactic Field Model 

From the sources in the simulation the probability vec- 
tors for extragalactic propagation are calculated for a user- 
defined discretization of the energies and directions. The 
calculation is separated into C++ classes for the prop- 
agation and the energy loss, each based on an abstract 
interface. The abstract interfaces are implemented as sub- 
types for the random- walk propagation in turbulent fields, 
respectively the described energy loss for proton and iron 
UHECR. This polymorphic design enables users to modify 
and extend the individual components independently. 

Galactic Field Model 

For an angular resolution of the discretization better 
than ^ 1° matrices of about 50000 x 50000 elements are 
needed. However, as in typical galactic field models par- 
ticles from most directions are not distributed over the 
whole sky, the matrices Ci are only sparsely populated. 
The lenses for the galactic fields are consequently imple- 
mented using a common linear algebra library which fea- 
tures sparse matrices [22 . This enables calculation of 
eq. 15 with reasonable consumption of resources. PAR- 



SEC includes tools for generation of the lenses from back- 
tracking data from the CRT program [9]. 

The galactic lenses are independent of the PARSEC 
module for the extragalactic propagation and can be used 
to calculate the deflection of individual cosmic rays. Spline 
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interpolation and numeric integration routines used in the 
program are taken from the GNU Scientific Library [23] . 

3. Results 

3.1. Energy spectrum 

The energy spectrum of cosmic rays is a key distribu- 
tion for comparisons of models with observations. From 
PARSEC simulations the energy spectra dN/dE(Ei) are 
calculated using dN/dE(Ei) = L • with normaliza- 

tion factor L representing the absolute scale to the relative 
luminosities of the sources in eq. [T] 
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Figure 1: Energy spectra of UHECR generated by PARSEC and ob- 
served spectra reported by the Pierre Auger Observatory [24] and 
the HiRes experiment |25j . The simulation results are obtained from 
sources with spectral indices 7 = —2.7 and 7 = —2.0 in an extra- 
galactic magnetic field with B = 1 nG and A = 1 Mpc. The shaded 
regions indicate the spread in the simulated realizations. 

In Figure [l] energy spectra obtained with the PARSEC 
program are compared with the observed energy spectra 
of the Pierre Auger Observatory [24] and HiRes experi- 
ment [25 . The spread of the energy spectra of 50 realiza- 
tions with different source positions and the mean energy 
spectrum are shown for two different spectral indices 7 of 
the injection spectra. The simulated spectra were gener- 
ated using isotropically distributed sources with a density 
of 10 -5 Mpc -3 and an extragalactic field of strength B = 1 
nG with correlation length A = 1 Mpc. The probability 
maps p i have been calculated for 100 log-linear spaced en- 
ergies Ei from 10 18,5 eV to 10 20,5 eV in the simulation. 
The normalization factor L has been set to match the re- 
sult from the Pierre Auger Collaboration at an energy of 
22 EeV. 



3.2. Particle horizons 

From the maximum injection energy E max and the 
model of the energy loss a particle horizon can be derived 
which corresponds to the maximum linear distance a par- 



ticle can originate from. Following eq. 12 this horizon de- 
pends on the extragalactic magnetic field model and the 
observed particle energy. In Figure [2] the horizon for differ- 
ent field strengths is shown as a function of the energy for 
sources with a maximum injection energy of E max = 1000 
EeV. 
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Figure 2: Horizon d for protons with observed energy E. The curves 
show d for several strengths of a turbulent magnetic field B with 
correlation length A = 1 Mpc. 

In case of a non-zero extragalactic magnetic field, the 
horizon is reduced as the trajectories are elongated. The 
effect is visible below 20 EeV for a 1 nG field and at higher 
energies for stronger fields. Without scaling the energy 
losses with red-shift z, the horizon is not reduced for lower 
energetic particles, as the energy loss only depends on the 
particles energy and not the time. The energy thus mono- 
tonically increases with the linear distance. 

For an energy of 100 EeV the distance from within 90% 
of the UHECR flux originates from is calculated to dgo = 
45 Mpc in case of zero magnetic field and a spectral index 
of the sources 7 = 2.0. This is « 30% lower than the result 
obtained from simulations with CRPropa [6] and a known 
feature of the continuous energy loss approximation [26] . 

3.3. Mean Deflection in Magnetic Fields 

The expectation value for an angle following a Fisher 
distribution with concentration parameter k is calculated 

as 
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Figure 3: Mean deflection of cosmic rays in (a) turbulent extragalactic fields and (b) BSS_S and JF galactic field model. 



with Io representing the modified Bessel function of order 
0. From this the mean deflection in the extragalactic field 
can be calculated for the parameters of the magnetic field 
and source model using eq. [6] 

In Figure [3] a) the mean deflection of protons from a 
source at 10 Mpc distance in extragalactic magnetic fields 
with three different strengths and a coherence length A = 1 
Mpc are shown as function of the energy. A mean deflec- 
tion of 90° corresponds to an isotropic arrival distribution 
of the UHECR. 

The mean deflections 6 in the galactic field can be 
directly calculated from the lens Ci as 

© oc ^2 ^2 lm > n arccos ( e ™ • e n) (18) 

n m 

with e n , em being the unit vector in the direction of pixel 
m, respectively n. In general the mean deflection in the 
galactic field depends on the source configuration due to 
the suppression of individual regions by the galactic lens. 
To get a mean deflection which is independent of the source 
configuration, we independently normalized every column 
n of Ci for the calculation of the mean deflection so that 
||(£i)n||i = 1- I n Figure [3] b) the mean deflection from 
lenses for the JF Model [27] and the BSS_S model [18] [28] 
of the galactic field are displayed. For the BSS_S model a 
field normalization of Bq = 0.48 fiG and scale heights of 
z\ = 0.95 kpc and Z2 = 4.0 kpc were chosen. 

3.4- Exemplary Probability Distributions 

To demonstrate the capabilities of PARSEC we gener- 
ated an exemplary result with source positions taken from 
the 12th edition of the catalogue of quasars and active 



galactic nuclei (AGN) by Veron-Cetty and Veron [29]. Ev- 
ery AGN of the catalogue up to a distance of 1000 Mpc 
has been considered. For the extragalactic field we chose a 
field strength B = 3 nG and a correlation length of A = 1 
Mpc. The resulting probability maps from the extragalac- 
tic propagation p l eg are shown in the top row of Figure 
for the two different energies £i = 5 EeV (Fig. [I] a) ana 
E 2 = 75 EeV (Fig.gb). 

The second row of Figure [4] (c,d) displays the same 
probability maps after application of the galactic lenses 
for a BSS_S model of the galactic field with a normaliza- 
tion of Bo = 0.48 fiG and scale heights of z\ = 0.95 kpc 
and Z2 = 4.0 kpc. The lenses have been created by back- 
tracking 10 6 protons with the CRT program [9] for 100 
log-linear spaced mono-energetic simulations from 10 18,5 
eV to IO 20,5 eV. The lenses and probability vectors have 
been discretized into 49,152 equal area pixels following the 
scheme implemented in the HEALpix software [30] . 

The third row of Figure [4] (e,f) shows the probability 
map after application of a purely geometric detector ac- 
ceptance of a typical earth-bound observatory [51] . 

3.5. Performance 

The exemplary simulation was performed in 6,690 sec 
using a single core of a Lenovo Thinkpad T400 notebook 
with 4 GB RAM and an Intel Core 2 Duo P8600 2.4 
GHz CPU. The notebook has been benchmarked with a 
SPECfp_base2006 rate of 12.2 and SPECint_base2006 rate 
of 15.1 [32]. Peak memory usage of the program was 0.5 
GB. The size of the BSS_S magnetic lens on disk is 262 
MB. 
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Figure 4: Hammer projected probability density maps for an observed energy of E = 5 EeV (left column) and E = 75 EeV (right column) in 
galactic coordinates. Black indicates higher, white lower probability on an arbitrary scale. From top to bottom the panels show: Probability 
density before applying the galactic lens (a,b), probability density after application of the galactic lens (c,d), probability density after application 
of a typical detector acceptance (e,f). 



The performance can be improved by approximating 
the individual contribution from distant sources as an i- 
sotropic component. This approximation introduces an 
error e to the extragalactic probability vector p 1 . An ex- 
ample realization with isotropic source distribution with 
10 -4 sources per Mpc 3 up to 1000 Mpc is used to demon- 
strate the effect. In Figure [5] the resulting error e = 
WPunbiased-Pbiasedh is displayed as a function of the com- 
putation time for various values of the cut-off parameter 



a (see Appendix A.2). Here a represents a measure of the 
strength of the anisotropic signal contribution not consid- 
ered in the calculation. By this method the computation 
time can be reduced by an order of magnitude for a = 0.05, 
which introduces an uncertainty less than 1%. 

4. Conclusions 

We have presented a parametrized Monte- Carlo gener- 
ator for fast creation of simulated UHECR datasets. The 



code includes effects of proton propagation through ex- 
tragalactic and galactic magnetic field models as well as 
energy losses during propagation. Also implemented is an 
extension to a simple model of iron propagation. Arbi- 
trary source models can be easily tested with the code. 
The generation of datasets is sufficiently fast to enable 
extensive parameter scans of the models. The modular 
design of PARSEC combined with the graphical steer- 
ing is suited for user extensions of the code and imple- 
mentation of user defined models. The source code of 
PARSEC as well as exemplary magnetic lenses are avail- 



able under a GNU General Public License from http: 
/ / www . physik . rwth-aachen . de/parsec , 
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which depends on the extragalactic probability density 
p eg1 or the individual configurations of the source and ex- 
tragalactic propagation models, respectively. 

If p eg is known, the uncertainty can be calculated, as 
we can approximate the true lens £ by the mean of indi- 
vidual realizations. In the following calculations only two 
realizations £\ and £2 of £ are used in order of clarity. 
For two realizations the true lens can be approxim ated 



IO" 1 

Computation time t(a) [t(a = 0)] 



Figure 5: Computation time for the exemplary scenario and resulting 
error e for different choices of the cut-off parameter a. Times are 
given in units of the time needed for the full calculation. 



Appendix A. Appendix 

Appendix A.l. Uncertainties of Galactic Lenses 

From two realizations of the same model for the galac- 
tic magnetic field an upper limit of the introduced error 
can be derived as follows. For eq. 15 it is ||p 5<J|i < 1 
as individual regions of the sky are suppressed. However, 



from £ G 



with Z m>n we can generate 
/ h,i ••• h,N \ 



£ = 



In,i 
si 



In,n 
s N 



(A.l) 



with £ G R^+ix^ and s n = 1- ||Z n ||i such that £p eg = 
p Q with p Q G R Ar+1 and p^ — (pi, • • -Pn-> s). s n represents 
the suppression the UHECR flux from the extragalactic 
direction n and s = ^ s n the total suppression of p eg by 
£. By this definitions it is ||p ||i = 1. 

Let C\ be a realization of the 'true' lens £, then ap- 
plication of £\ in eq. 15 introduces an uncertainty Sp 



A • Pea =Po + S P 



(A.2) 



as £ = \{£\ +£2)- Using this we substitute £\ in eq. 
which yields 



A.2 



S£p ea = Sp 



(A.3) 



£\ — £2 



with S£ 
For 

Sp = e • p Q (A.4) 

resembling a uniform uncertainty on the sky and unknown 
p eg we can estimate e by applying the unity norm ||.||i to 



eq. A.3 Using the Cauchy — Schwarz inequality this reads 



M < ^ll^lli 



(A.5) 
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as WPegh = \\Po\\l = L 

The definition of the unity norm reads ||<S£||i = max^ 
with Sl n being the n-th column vector of 6£. 



WSL 



Using 6£ as in eq. A.3 



and 



£ as in eq. 

= I 1 -I 2 



A.l 



SL 



(Sli in ,'"l Nin ,5s n ) With 5l m ,n = lL,n l m,n bein g tne 

difference of the elements of the matrices and Ss n = 
— being the difference of the corresponding suppres- 
sion factors. Thus we can write 



\Sl n \\i = 



I 1 



^m,n\ ~l~ I* 



(A.6) 
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(A.7) 

= ||5£ n ||i + |||/:i|| 1 -||£ 2 ||i| (A.8) 
using the definitions of the suppression factors. This yields 



e< J||«5i||i = Jmaxdl^ll! 

2 2 n 
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(A.9) 

as an upper limit of the uncertainty of the lens. 

The formalism can be extended to give an upper limit 
of the uncertainty in individual directions by substitution 



of the scalar e in eq. A.3 with a diagonal matrix S with 



elements e m5m being the relative uncertainty of the prob- 
ability in pixel n. Following the same calculation steps as 
above this yields 

1 



£ Po = -e m -5C-p e 



(A.10) 



with beeing the e m unit vector in direction m. Conse- 
quently this transforms to 



1 - 1 - T 

< " SC Wl = 2 m n X6lm 



(A.ll) 
(A.12) 



~ 1 

with row vector Sl m as upper limit of the uncertainty in a 
specific direction. 

From two realizations of the lenses used in section [3] we 
found a maximum uncertainty of 23% for a UHECR en- 
ergy of 10 18 " 5 eV with an typical uncertainty for individual 
pixels of about 2%. The maximum uncertainty above en- 
ergies E = 10 19 - 5 eV is less than 1%. 

Appendix A. 2. Performance Boost and Simulation Preci- 
sion 

The individual cosmic ray flux from many sources at 
large distances to the observer add up to give an almost 
isotropic contribution. The computation time spent to cal- 
culate this background can be eliminated by aborting the 
calculation for every individual pixel and adding the total 
isotropic background contribution to every pixel. By this 
we introduce an error e to p l eg . To check if the upcoming 
contributions are isotropic and decide whether to abort 
the detailed calculation we proceed as follows: first we di- 
vide the sources into 20 distance bins and calculate the 
contribution to p\ g from all sources in the first bin. We 
calculate a factor a = L upc -(mdixp l eg — m.mpl g )/(maxp l eg + 
minpg^) with L upc = ^2 fs being the integrated luminos- 
ity of all sources further away. If a is lower than a given 
cut off value, the upcoming luminosity is considered to be 
isotropic as the contribution from sources further away is 
more isotropic than from nearer sources. The flux from the 
upcoming bins is integrated and added once to all pixels 
in p l eg . If a exceeds the selected cut off value we proceed 
with the sources in the next bin. 
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